!! CONSTANT PARAMETERS

real(dp), parameter :: pi = atan(1.0_dp)*4.0e0_dp
real(dp), parameter :: Avogadro = 6.02214076e26_dp !![number/kmol]
real(dp), parameter :: Boltzmann = 1.380649e-23_dp !![J/K]
real(dp), parameter :: R0 = Avogadro*Boltzmann !! J/(kmol*K)
